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Abstract 

A one-dimensional reaction-diffusion model consisting of two species of 
particles and vacancies on a ring is introduced. The number of particles in 
one species is conserved while in the other species it can fluctuate because 
of creation and annihilation of particles. It has been shown that the model 
undergoes a continuous phase transition from a phase where the currents 
of different species of particles are equal to another phase in which they 
are different. The total density of particles and also their currents in each 
phase are calculated exactly. 

Recently one-dimensional reaction-diffusion systems have received much atten- 
tion because they show a variety of interesting critical phenomena such as out- 
of-equilibrium phase transitions A simple system of this type, which has 

been studied widely in related literatures, is the Asymmetric Simple Exclusion 
Process (ASEP) [5]. In this two-states model the particles are injected from 
the left site of an open discreet lattice of length L. They diffuse in the system 
and at the end of the lattice are extracted from the system. It is known that 
depending on the injection and the extraction rates the ASEP shows different 
boundary induced phase transitions. Non-equilibrium phase transition may also 
happen in the systems with non-conserving dynamics : 4, 5] . For instance in [5 
the authors investigate a three-states model consists of two species of particles 
besides vacancies on a lattice with ring geometry. The dynamics of this model 
consists of diffusion, creation and annihilation of both species of the particles. 
They have found that the phase diagram of the model highly depends on the 
annihilation rate of the particles. By changing the annihilation rate of the par- 
ticles, the system transfers from a maximal current phase to a fluid phase. The 
density of the vacancies changes discontinuously from one phase to the other 
phase. In present paper we introduce and study a reaction-diffusion model on 
a discrete lattice of length L with periodic boundary condition. Besides the va- 
cancies there are two different types of particles in the system. Throughout this 
paper the vacancies and the particles are denoted by E, A and B. The dynamics 
of the system is not conserving. The particles of type A and type B hop to the 
left and to the right respectively. The total number of particles of type B is a 
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conserved quantity and assumed to be equal to M. The density of these parti- 
cles is defined as ps = j^- In contrast, the total number of particles of type A 
is not a conserved quantity due to the creation and annihilation of them. Only 
the nearest neighbors interactions are allowed and the model evolves through 
the following processes 

Mb — > %A with rate 1 
$B — > 50 with rate 1 

AB — ► BA with rate 1 (1) 
A$ — ► 00 with rate us 
00 — ► A% with rate 1 

As can be seen the parameter us determines the annihilation rate for the particles 
of type A which besides the number of the particles of type B i.e. ps arc the free 
parameters of the model. One should note that the annihilation in our model 
only takes place for one species of particles. Our main aim in the present work 
is to study the phase diagram of the model in terms of us and the density of the 
B particles. 

In our model if one starts with a lattice without any vacancies the dynamics 
of the model prevents it from evolving into other configurations consisting of 
vacancies. In this case the system remains in its initial configuration and the 
steady state of the system is trivial. In order to study the non-trivial case 
we consider those configurations which have at least one vacancy. In order to 
find the stationary probability distribution function of the system we apply the 
Matrix Product Formalism (MPF) first introduced in [3] and then generalized in 
6 . According to this formalism the stationary probability for any configuration 
of a system with periodic boundary condition is proportional to the trace of 
product of non-commuting operators which satisfy a quadratic algebra. In our 
model we have three different states at each site of the lattice associated with 
the presence of vacancies, the A particles and the B particles. We assign three 
different operators E, A and B to each state. Now the unnormalized steady 
state probability of a configuration C is given by 

P(C) = ^-Tr[l[X i } (2) 

i—l 

in which Xj = E if the site i is empty, Xj = A if the site i is occupied by a 
particle of type A and Xj = B if it is occupied by a particle of type B. The 
normalization factor Zl in the denominator of ^ is called the partition function 
of the system and is given by the sum of unnormalized weights of all accessible 
configurations. By applying the MPF one finds the following quadratic algebra 
for our model 

AB = A + B 
AE = E 

EB = E 1 ' 

E 2 = wE. 

Now by defining E = V) (W| in which (W|V) = 1 one can simply find 

AB = A + B 

(W\B = (W\ 1 > 

E = us\V){W\. 
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The first three relations in (HJ) make a quadratic algebra which is well known 
in the related literatures. It is the quadratic algebra of the ASEP when the 
boundary rates are equal to ontJE This algebra has an infinite dimensional 
representation given by the following matrices and vectors 
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,B = A T ,\V) 
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V : J 



,(W\ = \V)' 



(5) 



in which T stands for transpose. In order to find the phase structure of the 
system one can calculate the generating function of the partition function of 
the system and study its singularities. In what follows we first calculate the 
grandcanonical partition function of the system Zl (£) by introducing a fugacity 
£ for particles of type B. We then fix the fugacity of the B particles using the 
following relation 

£ dlnZ L (Q 

The grandcanonical partition function of the system Z(£) can now be calculated 
from ^ and is given by 



(7) 



i=l 



One should note that the operator B in is replaced with the operator £B in 
which £ should be fixed using ©. As we mentioned above the stationary state of 
the system without vacancies is a trivial one, therefore in ([7J we have considered 
those configurations with at least one vacancy. The generating function for 
Zl(0 can now be calculated using ([!]). Using the same procedure introduced 
in [7] and after some straightforward algebra one finds 



L=0 



\ L - L z L (t) 



l-wtf&A) 



in which 



U(^X)^J2X L+1 (W\(A + ^B) L \V). 



(8) 



(9) 



L=0 



The convergence radius of the formal series |(HJ) which is the absolute value of 
its nearest singularity to the origin can be written as 



R(0 = lim Z L (0~ 

L — >oo 



(10) 



This is also the inverse of the largest eigenvalue of Zl (£) . In the large L limit 
using ([6]) this results in the following relation 

d . 1 



(11) 



1 The operator A and B should be regarded as the operators associated with the particles 
and vacancies respectively. 
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Therefore one should only find the singularities of ([5]) and decide in which region 
of the phase diagram which singularity is the smallest one. In order to find the 
singularities of (JSJ) one should first calculate £/(£, A). This can easily be done 
by noting that the matrix A + £B defined as 



A + £B 



/ i + € i 

€ i + e 

o £ 
o o 





1 

i + £ i 



\ 



(12) 



satisfy the following eigenvalue relation 

(A + £B)|0) = (1 + £ + 2v/£C OS (fl))|0) 
for — 7r < < 7r in which we have defined 



/ Sin(0) 
g-l 2 Sin{26) 

£Sm(36») 
C 3/2 S , m(46») 
£ 2 Sin{56) 

V 



(13) 



(14) 



/ 



Considering the fact that ^■Sin(9)Sin(n9) = <5i iM one can easily see that 

d0 



-Sin(0)\B). 



(15) 



Now using ([T3 ^) -([I5] ) and the fact that (W\8) — Sin{6) one can easily show that 
A) is given by 



U(£, A) 



d0 



ASm((9) 2 



(16) 



w tt l-A(l+e + 2^Co S (^)) 
. The integral l|16p can easily be calculated using 



which is valid for A < 

(i+V£) 

the Cauchy residue theorem by noting that it has three poles in the complex 
plane. Two of them are inside the contour which is a circle of unite radius 
around the origin and one is outside it. After some calculations one finds 



1 - A - A£ - y/(A + Ag - l) 2 - 4A 2 g 
2A£ ' 



(17) 



Now that [/(£, A) is calculated one can simply find the singularities of (jSJ. It 
turns out that (JHJ) has two different kinds of singularities: a simple root singu- 
larity Ri = n+ 1 ,jW< J j+?n which come from the denominator of (jHJ) and a square 
=-=— . Therefore the model has two different phases. 

The relation between the density of B particles and their fugacity in each phase 
should be obtained from (fTTI) . In terms of the density of the B particles ps we 



(i+w)(w+£) 
root singularity i? 2 = 



find Rt 



2 s - and i?2 = (1 — Pb) ■ Two different scenarios might happen: 
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defining ui c = 1 PB pB we find that for u> > uj c the nearest singularity to the origin 
is Ri and for uj < ui c it is i?2. The density of the vacancies can be calculated 
quit similar to that of the B particles. It is given by 

pE=uj iL in w) (i8) 

in which is again the nearest singularity to the origin in each phase. The 
density of the A particles is in turn pa = 1 — Pb~ Pe- Let us now investigate the 
current of the particles in each phase. Noting that the configurations without 
vacancies are inaccessible, the particle current for each species is obtained to be 

= Tr[(£AB + AE)(A + £B + E) £ ' 2 - (£AB)(A + gB)^- 2 ] 

Tr[(A + £B + E) L - (A + £B) L ] { J) 

= Tr[(£AB + £EB)(A + gB + E) £ ' 2 - (£AB)(A + gB) 1 - 2 ] 
B Tr[(A + eB + E) i -(A + ^B) i ] ' l j 

These relations can be simplified in the thermodynamic limit and one finds 

J A = R(0(l + {£-l)p A ) (21) 
Jb = R(m + (1-0PB). (22) 

As we mentioned above for uj < uj c the nearest singularity to the origin is 
always i?2- The particle currents in this case are equal and we find Ja = Jb = 
Pb(1 — Pb)- In contrast for uj > u c the nearest singularity to the origin is always 
Ri and it turns out that the currents are not equal. We find J a = [j^p an d 
Jb = Pb{^~Pb)- I n what follows we bring the summery of the results concerning 
the phase structure of the system 

!pA = l-pB Ja= Pb(1- Pb) 
Pb = Pb Jb = Pb(1-Pb) 
Pe = 

and 

!PA = I+ZJ J A = (1+^)2 
Pb = Pb Jb = Pb(1 - Pb) ■ 

PE = T^-PB 

As can be seen for to < co c the density of vacancies is equal to zero which means 
there are only A and B particles on the lattice. Since the density of the B 
particles is fixed and equal to ps the density of A particles should be 1 — ps- In 
this case, according to (UJ , both A and B particles have simply ASEP dynamics 
and therefore their currents should be of the form p(l — p) and that J a = Jb- 
This is in quite agreement with our calculations for J a and Jb- On the other 
hand for u> > ui c the density of vacancies on the lattice is no longer zero. At the 
transition point uj c the density of the vacancies is zero but it increases linearly 
in this phase. In terms of the density of the vacancies the phase transition is a 
continuous transition. In this phase for u> < 1 we always have J a > Jb while 
for uj > 1 we have J a > Jb for pb < j^j and J a < Jb for < pb < j^- 
In this paper we studied a three-states model consists of A and B particles 
besides the vacancies. The A particles are created and annihilated which is 
controlled by uj while the B particles only diffuse on the lattice and have a fixed 
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density ps- We found that the system have two phases depending on ps and u>. 
The current of the B particles is always a constant — pb) throughout the 
phase diagram while for the A particles it is given by different expressions in 
each phase. As a generalization one could also consider a more general process 

A$ — ► A with rate a 
$B — ► B$ with rate 

AB — > BA with rate 1 (23) 
A$ — ► 00 with rate A 
00 — > A$ with rate A' 



with the quadratic algebra given by 

= A + B 

r\ - \ v\ 

(24) 



AB = A + B 

a A\V) = \V) 
f3 (W\B = (W\ 
E=Tfe\V)(W\. 



and apply the same approach used in present paper to study its phase diagram. 
This is under our investigations. 
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